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FOREWORD 

This  technical  summary  report  was  prepared  by  the  Research  Department  of  Allison  Division 
of  General  Motors.  It  describes  the  work  accomplished  under  Contract  Nonr-4 104(00)  for  the 
period  1  May — 31  August  1065. 

The  program  was  sponsored  by  the  Office  of  Naval  Research  under  the  direction  of  Mr.  J.  A. 
Satkowski. 

The  report  was  written  bv  Dr.  Richard  T.  Schneider  and  Dr.  Horst  E.  Wilhelm.  Other  mem¬ 
bers  of  the  Allison  Research  team  who  contributed  to  the  work  were  R.  L.  Koch,  D.  L.  Clingman, 
R.O.  Whitaker,  R.  L.  Holtman,  andJ.A.  Eriksen,  Jr. 

Heater  design  and  fabrication  supervision  was  provided  by  Prof.  W,  Kluge  and  Dr.  W.  Bloss 
of  the  Technical  University,  Stuttgart,  Germany,  consultants  to  Allison  Division  of  General 
Motors. 
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I.  INTRODUCTION 

This  report  describes  the  progress  on  Contract  Nonr-  .  104(00)  for  the  period  1  May  —  31  August 
1965.  Design  and  operation  of  a  relativelv  large,  closed-loop  MFD  power  generation  research 
device  are  discussed. 


The  objective  of  the  effort  during  the  report  period  was  to  collect  information  which  would  per¬ 
mit  determining  the  feasibility  of  closed-loop  MPD  power  generation  at  low  temperatures  —  i.e., 
1500°K.  It  was  felt  that  it  would  be  necessary  to  employ  a  rue  closed-loop  device  rather  than 
a  blowdown  device  since  the  blowdown  device  avoids  many  of  the  difficulties  necessary  for  study 
and  gives  data  of  a  different  type.  The  data  which  would  be  obtained  by  the  blowdown  device 
would  require  too  much  extrapolation.  Also,  with  the  blowdown  device  ine  observation  time  is 
short  and  the  electrodes  are  cold  rather  than  at  operational  temperature. 

The  difficulties  encountered  as  n  result  of  (1)  investigat'ons  described  herein,  (2)  previous 
investigations,*  ’  (3)  investigations  emploving  a  smaller  closed -loop  device,  and  (4)  investiga- 
tions  using  three  other  existing  true  closed-loop  devices  -  -  were  of  a  nature  which  could  not 
be  detected  in  a  blowdown  device. 

The  main  difficulties  encountered  in  the  closed-loop  device  are  summarized  as  follows; 

•  Corrosion  of  the  components  caused  by  cesium  and  by  traces  of  oxygen  and  nitrogen 

•  Obtaining  a  smooth,  steady  injection  of  helium  over  extended  time  periods 

•  Separation  of  helium -cesium 

•  Partial  or  complete  electrical  shorting,  caused  by  conducting  deposits,  between  electrodes 
and  between  electrodes  and  ground 

•  Design  of  a  reliable  heater  capable  of  heatir.g  a  large  mass  flow  of  helium  without  intro¬ 
ducing  impurities 

Some  of  these  difficulties  have  been  overcome  and  some  need  further  work,  as  discussed  in 
detail  herein. 


^Superscripts  refer  to  references  in  Section  VIII. 
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GENERATOR  DESIGN 


DUCT 

Figure  1  shows  a  cross  section  of  the  new  M PD  generator  test  section.  The  walls  of  the  duct 
(center  portion  of  Figure  1)  are  formed  by  two  layers  of  ceramic  backed  with  a  layer  of  tanta¬ 
lum.  The  layers  of  ceramic  are  kept  from  moving  parallel  to  the  tantalum  wall  bv  stainless 
steel  brackets  in  the  cold  parts  of  the  duct.  Ceramic  spacers  ai  e  used  to  keep  ceramic  layers 
and  the  tantalum  wall  from  moving  toward  each  other.  The  two  ceramic  layers  are  kept  in 
place  by  the  shoulders  on  both  sides  of  the  ceramic  spacers.  These  spacers  form  two  cavities 
with  ceramic  shields  between  the  cavities.  Oversized  holes  in  the  shields  prevent  the  ceramic 
material  from  touching  the  electrode  stem  before  final  contact  is  made  in  the  colder  portion  of 
the  duct  where  the  electrode  stem  is  held  mechanically.  Mechanical  contact  between  the  elec¬ 
trode  and  the  ceramic  is  made  by  an  electrode  holder.  The  electrode  holder  is  sufficiently 
removed  from  the  hot  zones  of  the  duct  to  force  leakage  current  in  the  ceramic  to  go  through 
cold  ceramic  parts,  where  the  electrical  resistance  is  high. 

With  the  cavity  design  as  described,  it  was  possible  to  overcome  the  electrode  shorting  prob¬ 
lem  caused  by  deposits  of  conducting  material.  The  shield  plate  in  front  of  the  first  cavity 
prevents  most  of  the  material  from  being  deposited  inside  the  first  cavity.  To  short  the  elec¬ 
trode,  the  material  being  deposited  must  pass  over  two  more  shield  plates  and  two  more 
cavities  —  i.e,,  make  a  continuous  plating  —  before  an  electrical  contact  is  made.  The  photo¬ 
graph  of  the  ceramic  parts  in  Figure  2  illustrates  more  clearly  how  the  cavity  dosi~n  works. 
Figure  2  also  shows  how  the  inner  plates  are  held  in  place  and  how  the  outer  plates  arc  held 
apart. 

The  ceramic  parts  of  the  test  section  are  made  of  alumina.  In  recent  runs,  however,  the  inner 
plates  were  made  of  boron  nitride.  The  physical  size  of  the  channel  shown  in  Figure  2  is  as 
follows: 

•  Electrode  distance  —  7  cm 

•  Channel  width  (between  inner  plates)  —  2  cm 

>  Total  channel  length  (consisting  of  several  elements  as  shown  in  Figure  2)  — 50  cm 

A  cooling  channel  is  provided  in  the  generator  test  section  between  the  magnet  pole  pieces  and 
the  tan  Turn  wall  (re  Figure  1).  Cold  helium  gas  is  blown  into  this  channel  to  protect  the  mag¬ 
net  pole  pieces  from  heat.  This  gas  is  later  mixed  with  the  main  gas  stream.  The  magnet  pole 
pieces,  which  are  water  cooled,  form  a  part  of  the  outer  wall  of  the  duct  in  such  a  way  that 
they  provide  mechanical  strength  without  having  a  sealing  effect. 
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Figure 


1.  Cross  section  of  the  MPD  generator  test  section  (reference  EX-71131). 
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Figure  2.  Ceramic  parts  for  cavity  design.  Note  maimer  of  holding  inner 
plates  in  place  and  outer  plates  apart. 

The  electrode  distance  can  be  varied  by  interchanging  the  inner  ceramic  plates  of  the  duct. 

The  height  of  these  plates  determines  the  electrode  spacing. 

Figure  3  shows  the  side  view  of  the  generator  duct.  A  diagnostic  box  is  attached  to  each  end 
of  the  channel.  Each  box  has  a  straight  port  for  spectroscopic  measurements  and  an  inclined 
port  for  measurement,  of  me  electrode  surface  temperature.  For  additional  observations,  a 
pair  of  windows  is  mounted  directly  to  the  duct.  Optical  access  was  obtained  by  removing  one 
electrode  pair  and  using  the  stem  holes  as  apertures.  The  number  of  electrode  pairs,  which 
varies  between  57  and  67,  depends  on  the  electrode  spacing  used. 

In  the  present  buildup,  different  electrode  materials  arc  used  alternately.  These  art  pure 
tungsten.  Philips  Type  A,  and  Philips  Type  B.  The  segmentation  is  as  follows: 

•  Distance  between  center  Enes  of  neighboring  electrodes— 0.75  cm 

•  Electrode  width  —  0.5  cm 

•  Free  space  between  two  neighboring  electrodes  —  0.25  cm 
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The  exit  section  of  the  duct  is  identical  with  the  entrance— a  diagnostic  box  is  provided  and  may 
be  used  to  obtain  other  diagnostic  measurements. 

Figure  4  is  a  photograph  of  the  generator  duct.  The  magnet  pole  pieces  are  attached.  The 
observation  windows  are  also  visible. 

HEATER 

Tt  e  experience  gained  from  operation  of  the  small  heaters1  was  invaluable  in  designing  the  new 
heater,  as  described  in  the  following  paragraphs.  First,  it  was  impossible  to  support  a  heater 
element  with  an  insulator  when  the  heater  was  in  the  horizontal  position — heater  element  tem¬ 
perature  being  necessarily  too  high  for  existing  insulator  materials  to  withstand.  Consequently, 
the  heater  element  was  positioned  vertically  so  that  the  insulator  would  be  near  the  cooler  por¬ 
tion  only.  (The  heater  element,  when  positioned  vertically,  does  not  require  support  in  those 
areas  which  run  hot.) 

Secondly,  it  was  necessary  to  design  the  heater  element  so  that  very  little  machining  would  be 
required  in  going  from  the  raw  material  shape  to  the  final  shape.  Each  forming  process  in¬ 
vites  cracks  which  ultimately  cause  failure  during  operation.  Another  requirement  was  that 
the  heater  element  be  designed  to  keep  contamination,  in  form  of  metal  vapor,  in  the  gis  to 
a  minimum.  As  a  final  requirement,  it  was  necessary  that  the  heater  be  capable  of  shutdown 
from  maximum  to  zero  power  within  a  few  seconds. 

The  final  design  of  the  heater  which  meets  the  requirements  stated  previously  is  shown  in 
Figure  5.  Since  the  heater  was  designed  for  vertical  operation,  the  generator  also  operates 
in  the  vertical  position. 

As  shown  in  Figure  5,  helium  enters  through  the  helium  inlet  flange  at  room  temperature.  It 
flows  through  the  tungsten  radiation  shields  and  cools  them.  On  return  through  the  heater,  the 
helium  is  preheated  as  it  passes  along  the  tungsten  box  containing  the  heating  elements.  (The 
heating  elements  are  not  in  direct  contact  with  the  working  fluid.)  After  passing  the  tung.  ‘en  box, 
the  helium  enters  the  main  heater  where  it  is  brought  up  to  2000°K  by  radiation  heating. 

The  difference  between  gas  temperature  and  wall  temperature  in  the  heater  section  is  not  ex¬ 
pected  to  ,e  more  than  200°K.  All  high  temperature  parts  of  the  heater  are  made  of  tungsten 
to  minimize  metal  vapor  pressure. 

Following  are  the  technical  data  concerning  the  helium  heater: 

•  Maximum  oower  input  — 180  kw 

•  Design  current — 5001  amp/pbooe,  three-phase,  ac 

•  Design  voltage  —  20  v,  three-phase,  ac 
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Figure  3.  Side  view  of  the  MPD  test 
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Figure  3.  Side  view  of  the  MPD  test  section. 
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A  three-phase  regulated  transformer  is  used  having  the  primary  winding  for  440  v.  The  voltage 
output  on  the  secondary  is  0  to  3  5  v.  The  transformer  is  rated  at  2u(J  kva. 

Figure  0  shows  a  photograph  of  the  completely  assembled  MPD  power  generator  research  de¬ 
vice.  Design  data  for  the  MPD  power  generator  system  ai'e  presented  in  Table  1. 


Figure  4.  Generator  duct  with  magnet  pole  pieces  attached. 
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Table  1. 

MPD  power  generator  design  data. 


Working  fluid 
Pressure  range 
Seed 

Seeding  methods 
Mass  flow 

Heater  power  input — gross 

Heater  voltage 

Heater  current 

Duct  length 

Electrode  distance 

Channel  width 

Number  of  electrode  pairs 

Segmentation  of  electrodes 

Electrode  width 

Electrode  materials 


Magnetic  field  strength 


Helium 

0.5  to  2.5  atm 

Cesium 

Cesium  boiler  and  superheater 

Up  to  gm/sec  helium  with  two  pumps  parallel 

Up  to  15  gm/sec  helium  with  one  pump 

180  to  200  kw  (maximum) 

20  v,  three-phase,  ac 
3000  amp/ phase 
50  cm 
7.5  cm 
2  cm 
Up  to  67 

0.75  cm  (center-to-center) 

0.5  cm 

Porous  tungsten 
Philips  electrodes  type  A 
Philips  electrodes  type  B 
2.15  Tesla  (maximum) 


Cs -Supply 
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III.  M PD  GENERATOR  CHECKOUT  RUNS 

The  objective  of  the  checkout  runs  is  to  establish  calibration  curves  which  allow  predetermina¬ 
tion  of  the  operating  conditions  for  a  1  actual  cesium  run.  Also,  they  serve  as  an  endurance  test 
with  regard  to  determining  system  limitations  and  weaknesses,  thus  reducing  the  chance  for  a 
failure  during  an  actual  cesium  run. 

The  accumulated  running  time  of  the  new  system  is  228  hr  to  date.  Only  about  5%  of  this  time, 
however,  was  actual  cesium  runs.  The  calibration  curves  established  during  the  checkout  runs 
consist  of  heater  power  input  versus  gas  temperature,  with  mass  fbw  and  pressure  as  param¬ 
eters.  Thermocouples  and  optical  pyrometry  were  used  for  measurement  of  the  gas  tempera¬ 
ture. 

Most  of  the  checkout  runs  were  made  at  the  1250,  1340,  and  1450°K  levels.  The  longest  run 
made  at  a  specific  temperature  (1340°K)  without  changing  flow  rate  or  any  other  condition 
lasted  2  hr.  (Endurance  runs  at  lower  temperatures  were  often  made  overnight  without  obser¬ 
vation.) 

The  generator  duct  was  disassembled  only  once  during  the  accumulated  running  time  of  228 
hr — this  was  to  make  minor  modifications,  not  because  of  a  failure.  The  heater  has  not  been 
disassembled  at  all  to  date.  These  endurance  runs  show,  therefore,  that  there  are  no  materials 
and  design  problems  at  the  temperature  levels  encountered  to  date  and  indicate  that  temperature 
levels  of  1j00°K  are  possible. 

The  magnet  was  operated  up  to  a  maximum  field  of  21,  500  gauss.  This  value  was  measured 
with  the  test  section  in  place  employing  three  different  gauss  meters.  Two  of  the  gauss  meters 
had  solid-state  sensors  (Hall  generators)  and  the  third  had  a  rotating  coil.  The  values  obtained 
agreed  within  5%.  Since  it  is  impossible  to  have  a  gauss  meter  installed  when  running  under 
hot  conditions,  a  calibration  curve  of  magnet  current  versus  field  strength  was  established. 

This  curve  is  fed  into  an  analog  computer,  and  magnet  current  measured  during  a  run  is  con¬ 
verted  continuously  into  magnetic  field.  Therefore,  all  machine-plotted  curves  presented  in 
this  report  have  a  linear  scale  for  magnetic  field. 

The  effectiveness  of  the  cavity  design — designed  to  reduce  electrical  shorting  —  was  also  tested 
during  the  endurance  runs.  For  unseeded  runs  the  resistance  between  opposite  electrodes  was 
never  less  than  0.1  megohm.  This  was  also  true  for  most  of  the  seeded  runs.  Following  quick 
shutdowns  aftei  seeded  runs,  electrode  resistances  as  low  as  1000  ohms  were  observed.  Upon 
reheating,  however,  the  electrode  resistance  returned  to  normal  values.  It  is  concluded,  there¬ 
fore,  that  cesium  vapor  condensed  in  the  cavities  and  formed  a  conducting  film.  When  the  sys¬ 
tem  was  reheated,  this  film  evaporated.  It  was  possible  to  keep  all  other  metal  vapors  from 
forming  a  continuous  film. 
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During  the  unseeded  endurance  runs,  outputs  in  the  millivolt  range  were  observed.  These 
outputs  are  not  considered  to  be  MPD  generator  outputs  because  they  were  independent  of 
velocity  and  the  absolute  value  of  the  magnetic  field.  This  output  was  5  mv  per  electrode 
pair  at  room  temperature  and  went  up  to  50  mv  at  high  temperatures.  The  internal  impedance 
of  the  system  was  about  1  megohm;  therefore,  the  millivolt  outputs  never  disturbed  any  MPD 
measurements.  They  were  useful,  however,  for  instrumentation  checkout.  An  exact  explana¬ 
tion  for  this  millivolt  output  cannjt  be  given  as  yet.  It  can  be  stated,  however,  that  it  is  not 
an  electrical  pickup  from  the  heater  or  any  other  external  voltage  source,  since  it  has  also 
been  observed  with  all  power  shutoff. 

The  impurity  level  in  the  helium  gas  was  continuously  monitored  by  an  oxygen  analyzer  and  an 
electronic  hygrometer.  The  oxygen  level  was  maintained  near  1  ppm  and  the  water  vapor  nei.r 
25  ppm. 

In  addition  to  the  work  previously  described,  effort  was  initiated  to  measure  all  impurities 
with  a  mass  spectrometer.  The  gas  inlet  to  the  mass  spectrometer  was  designed  to  permit 
continuous  sampling  of  the  working  fluid.  Since  the  relationship  of  the  compu,  ition  of  the  gas 
in  the  mass  spectrometer  and  the  sample  to  be  analyzed  must  be  known,  a  viscous  leak  design 
following  that  of  Halsted  and  Niers  was  selected.  With  this  type  of  inlet  design,  the  sample  is 
introduced  from  a  high  pressure  region  (approximately  1  atm  or  more)  directly  into  the  mass 

_5 

spectrometer,  which  must  operate  at  approximately  10  torr  or  less. 

The  gas  inlet  consisted  of  a  type  316  stainless  steel  capillary  tube  1.5  meters  long  with  an  ID 
of  1.78  X  I0“4  meter.  A  constriction  was  placed  in  the  capillary  at  the  entrance  to  the  mass 
spectrometer  and  was  adjusted  so  that  the  pressure  in  the  mass  spectrometer  was  approximately 
1  X  IQ'5  torr.  A  resistance  heater,  consisting  of  thermocouple  wire,  w'as  placed  around  the 
capillary  to  prevent  condensation  of  the  working  fluid  constituents  within  the  capillary  tube. 

The  mass  spectrometer  was  a  General  Electric  Model  22  PC  Partial  Pressure  Analyzer  with 
a  Model  22PT  110  type  304 L  stainless  steel  analyzer  tube.  The  tube  was  modified  to  permit 
introduction  of  the  gas  sample  in  the  immediate  vicinity  of  the  ion  cage.  This  modification 
eliminated  the  uncertainty  resulting  from  the  long  path  diffusion  normally  used  in  this  tube. 

Halsted  and  Nier*  have  shown  that  a  multicomponent  gas  flowing  through  a  viscous  leak  of  the 
type  previously  described  follows  approximately  the  relationship 


lb  ‘  LMbj  Pb 


as  the  upstream  pressure  becomes  large. 
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In  the  equation, 

I  mass  spectrometer  current 

M  =  molecular  weight  of  component 

P  component  partial  pressure  in  main  gas  stream 

a  and  b  =  identify  components 

Systematic  surveys  of  gas  samples  were  made  covering  the  mass  number  range  from  2  to  132. 
A  cross-check  of  the  results  was  made  by  comparing  mass  number  32  data  with  information 
obtained  from  a  Lockwood  and  McLorie  Oxygen  Analyzer.  Ratios  of  readings  taken  at  different 
times  from  the  two  instruments  never  varied  by  more  than  a  factor  of  3. 

Recordings  showing  the  helium  and  oxygen  ion  currents  for  both  gettered  and  nongettered  sam¬ 
ples  were  made.  Partial  pressures  for  oxygen  and  helium  as  calculated  from  the  Halsted-Nier 
equation  are  as  follows: 

Gettered  Nongettered 

883  torr 
1.02  X  10"  1  torr 

1.15  X  10"4 


’lie  = 

917  torr 

PHe  = 

02  = 

6.5  X  10-^  torr 

po2  r 

>02 

_CL 

P°2 

— —  - 

7.08  X  10  0 

He 

PHe 
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IV.  HELIUM -CESIUM  RUNS 

He-Cs  runs  were  made  at  1250.  1340,  and  1450°K.  The  He  mass  flows  varied  between  10  and 
13  grn/sec.  The  pressure  in  the  duct  was  varied  between  1  and  1.3  atm.  The  r.aximum  value 
of  the  magnetic  field  was  2.  15  Tesla. 

Figure  7  shows  a  plot  of  open  circuit  voltage  versus  magnetic  field — the  theoreti  al  vbd  (voltage, 
magnetic  field,  electrode  distance)  curve  is  indicated.  The  leakage  resistance  for  the  run 
plotted  in  Figure  7  was  1  megohm.  Also,  the  run  was  made  with  an  extremely  high  seeding 
ratio  C>2'7o).  However,  as  soon  as  the  seeding  ratio  was  reduced  to  a  moderate  level  and  the 
leakage  resistance  dropped  to  0.  1  megohm,  the  saturation  effect  observed  in  earlier  runs  re¬ 
occurs.  To  date  it  has  not  been  possible  to  distingmsh  which  of  the  two  effects  (leakage  resis¬ 
tance  and  seeding  ratio)  causes  the  saturation  effect.  Additional  runs  are  required  before  a 
decision  can  be  made.  Figure  8  shows  a  typical  curve  of  voltage  versus  magnetic  field. 

These  curves  are  reproducible,  as  has  been  shown  b'  numerous  runs. 


Figure  7.  Oper.  circuit  voltage  vs  magnetic  field— seeding  >2°'. 
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Figure  8.  Reduced  seeding  ratio  and  100  kilohm  leakage  resistance. 

In  the  He-Cs  runs,  the  open  circuit  voltage  has  been  reduced  considerably.  The  magnet n  field 
has  a  saturation  effect,  and  as  the  magnetic  field  increases  above  approximately  1.5  Fes  la,  the 
open  circuit  voltage  decreases.  Figure  8  was  plotted  directly  with  an  x-y  plotter  during  the  run. 

Figure  9  was  plotted  to  show  the  pronounced  decrease  in  open  circuit  voltage  in  the  2-Tesla 
region. 

To  show  the  saturation  effect  with  respect  to  the  influence  of  the  Hall  effect,  open  circuit  voltage 
was  recorded  using  different  electrode  lengths.  The  electrode  length  was  increased  by  connect¬ 
ing  cle<  trode  segments  together  electrically.  The  results  of  this  run  are  shown  in  Figure  10. 

I'lie  tes*  starts  with  curve  No.  3,  which  indicates  that  electrode  No.  3  was  connected  to  the  x-y 
recorder.  File  next  curve,  called  No.  3+4,  indicates  that  electrodes  No.  3  and  4  ate  connected, 
forming  one  electrode  twice  the  length  of  the  one  shown  in  curve  No.  3,  etc.  The  scale  fat  tot 
was  ch  tnged  tiv  a  factor  of  two  for  the  lower  three  curves  sn  Figure  lft. 
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As  shown  in  Figure  10,  lengthening  the  olcctrc  •o  altered  the  shape  of  the  curve  only  slightly. 
This  slight  change  might  be  attributed  to  the  varying  load.  Therefore,  it  can  be  concluded  that 
the  Hall  effect  is  not  the  dominant  cause  of  the  saturation  effect.  Since  each  electrode  pair 
has  a  different  leakage  resistance,  and  since  these  leakage  resistances  became  parallel 
when  the  electrodes  were  lengthened,  each  curve  was  taken  at  a  different  effective  load.  Tin 
effect  of  different  loads  on  the  saturation  curve  is  shown  in  Figure  8. 

As  mentioned  previously,  the  saturation  effect  can  be  overcome  bv  increasing  the  seeding  ratio. 
Figure  11  shows  two  curves  —  one  taken  with  a  seeding  ratio  <2%  and  one  taken  with  a  seeding 
ratio  >2%.  The  curve  with  the  smaller  seeding  ratio  shows  saturation,  while  the  curve  with 
the  larger  seeding  ratio  shows  no  saturation,  although  it  is  not  considered  to  be  a  good  curve 
because  the  cesium  mass  flow  was  fluctuating.  Thus,  the  saturation  effect  shows  up  only  on 
runs  with  moderate  or  small  seeding  ratios. 


Seed  >  2% 


Seed<  2% 


Tesla 
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Figure  11.  Saturation  overcome  by  increasing  the  seeding  ratio. 
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Vo  date,  it  lias  not  been  possible  to  run  l-Y  curves  in  tne  nonsaturated  operating  mode  ( Figure  7). 
Howevei  ,  a  number  of  1-V  curves  were  obtained  in  the  saturated  mode  (Figure  8). 

A  typical  I-V  curve  for  the  saturated  mode  is  shown  in  Figure  12.  The  straightness  of  the  lines 
indicates  that  no  electron  heating  was  observed.  It  is  believed  that  in  the  nonsaturated  mode, 
however,  electron  heating  could  be  observed.  This  information  will  be  obtained  during  the  next 
series  of  tests. 

In  addition  to  the  x-y  plots,  the  data  from  the  runs  made  during  the  report  period  were  also 
recorded  with  an  eight-channel  strip  chart  recorder.  One  channel  was  used  for  magnetic  field, 
one  for  helium  mass  flow,  and  the  remaining  six  were  used  to  plot  the  time  behavior  of  six 
different  electrode  pairs  selected  from  all  parts  of  the  generator  duct.  Figure  13  shows  a 
typical  strip  chart  recording.  The  saturation  effc  t  is  shown  by  the  coincident  lines  of  maxi¬ 
mum  magnetic  field  and  minimum  voltage  output  of  the  electrodes.  Equally  important  in  the 
recording  is  the  fact  that  the  shape  of  the  saturation  curve  does  not  change  when  going  from 
one  location  in  the  channel  to  another.  This  also  indicates  that  no  motion- induced  nonequili¬ 
brium  ionization  was  observed. 


Figure  12.  Typical  I-V  curves  in  saturated  mode. 
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The  electrodes  for  the  runs  were  numbered  progressively  starting  with  No.  1  at  the  channel 
entrance.  In  most  of  the  runs,  output  voltage  increased  with  electrode  number,  as  shown  in 
Figure  14.  The  increase  in  output  voltage  is  a  result  of  an  increase  in  gas  velocity  due  to  vis¬ 
cous  effects.  A  theoretical  curve  which  takes  this  viscous  effect  into  account  is  also  shown  in 
Figure  14. 

Finally,  an  attempt  was  made  to  measure  Hall  voltages.  The  results  are  shown  in  Figure  15. 
The  voltage  between  neighboring  electrodes  was  of  the  same  magnitude  as  the  voltage  between 
opposite  electrodes  —  it  was  much  too  low  to  cause  a  breakdown  between  neighboring  electrodes. 
When  connecting  the  electrodes  in  the  Hall  mode,  the  Hall  voltage  observed  was  less  than  1  mv 
between  two  neighboring  electrodes. 
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Figure  14.  Comparison  of  theoretical  and  observed  curves  of  output  vs  electrode  numbers. 
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Figure  15.  8-channel  measurement  of  Hall  voltages 
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V.  DISCUSSION  OF  RESULTS 

Based  on  the  experimental  results  given  in  the  preceding  sections,  two  questions  remain  un¬ 
answered.  First,  why  is  the  theoretical  open  circuit  voltage  not  observed  in  the  two  operational 
modes  1  Second,  how  can  the  saturation  effect  he  explained?  This  section  discusses  these 
questions. 

It  lias  been  determined  experimentally  that  the  Hall  effect  does  not  cause  the  saturation  effect 
{Figure  10).  Other  mechanisms  that  could  cause  this  effect  must,  therefore,  be  found.  One 
of  the  mechanisms  could  be  a  velocity  profile  over  the  cross  section  of  the  channel,  as  shown 
in  Figure  16.  As  a  result,  the  generated  emf  would  vary  across  the  channel.  For  this  mech¬ 
anism  an  electric  field,  S',  is  impressed  on  the  channel  by  the  load  resist*,  or  th«  open  cir¬ 
cuit.  In  the  case  of  t*  ue  open  circuit,  E'  is  equal  to  the  peak  value  of  emf  in  the  center  of  the 
channel.  Since  the  electrode  extends  across  the  channel,  F'  is  also  extended  across  the  chan¬ 
nel.  However,  a  smaller  emf  counteracts  this  E'  in  the  proximity  of  the  channel  wall.  Thus, 
there  is  current  flowing  in  the  proximity  of  the  wall  driven  by  emf  onifwall  ~  Since 
emfwail  <  E'»  the  current  flows  in  the  opposite  direction  of  any  current  generated  in  the  center 
of  the  channel.  Figure  16  illustrates  this  condition. 


Electrode 


Figure  16.  Velocity  and  emf  profile 
across  channel. 
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For  the  ease  of  external  open  circuit,  a  current  as  previously  described  flows  through  an  in¬ 
ternal  resistance  Hjs  (profile).  This  resistance  can  be  of  the  same  order  of  magnitude  us  the 
internal  resistance  rtj  since  it  involves  the  same  gas  and  comparable  conductivity.  Figure  17 
shows  that  if  Rj  Rg,  the  open  circuit  voltage,  Vnb,scrved.  would  be  only  50%  of  the  true  open 
circuit  voltage  (which  is  approximated  in  Figure  7).  Under  these  conditions,  as  well  as  at 
very  low  magnetic  fields,  the  true  open  circuit  voltage  can  never  be  observed.  This  has  been 
confirmed  b’  laboratory  observation. 

Another  possibility  exists  for  internal  leakage.  From  Figure  1G,  consider  the  walls  to  be 
coated  with  a  layer  of  conducting  material  and  connected  to  one  of  the  electrodes.  The  resis¬ 
tance  between  electrodes  would  still  be  high.  Fvcn  though  the  impressed  electric  field,  F', 
now  lies  between  the  one  electrode  and  the  conducting  wall,  which  is  connected  to  the  other 
electrode,  the  resistance  Rjg  (wall)  of  this  gap  is  on  the  order  of  Rj  since  the  gas  is  the  same 
and  the  conductivity  is  comparable. 

Referring  to  Figure  18,  the  situation  is  shown  where  a  load  resistor  is  applied.  Since  gaseous 
resistances  are  rarely  constant  over  a  wide  voltage -cur rent  range,  it  is  very  unlikely  that 
Rjg  (profile)  and  Rjg  (wall)  will  remain  constant. 


r-nif 

Figure  17.  Schematic  for  internal 
shorting— open  circuit. 

t 


VoliSHrv(’d  '''rout 


4400-22 


28 


Allison 


Tb  question  concerning  whether  or  not  resistance  is  a  function  of  current  or  ernf  can  ho  an¬ 
swered  bv  tlie  1-V  curves  in  Figure  12.  The  curves  approximate  siraignt  lines,  indicating  that 
current  alone  will  not  change  Rp  RJS  (profile),  and  RJS  (wall).  Therefore,  the  change  in  R, 
and  Rjg  which  causes  the  saturation  effect  must  be  a  function  of  emf. 


Figure  18  (right  side)  shows  how  Rj^  (profile)  and  Rj^  (wall)  are  affected  hy  emf.  The  observed 
voltage  is 

RL 

^observed  cm*eff  ^  p 


As  previously  stated,  Rj  =  f  (emf).  Therefore,  it  may  be  assumed  that 
Rt  a  Eb  f  c 

Vff 


and 


ernfyfj‘  -  d  *  mf 
where  a,  b,  c,  and  d  are  constants. 


I  4400-23 

figure  18.  Schematic  for  internal  shorting  with  load  applied. 
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The  observed  data  may  be  matched  by  proper  selection  of  the  constants  a,  b,  c,  and  d.  The 
exact  values  for  >se  constants  will  be  established  and  reported  later. 

It  is  desirable  to  distinguish  between  the  effect  of  emf  on  (wall)  and  Rjg  (profile).  This  is 
possible  using  the  diagram  on  the  left  side  of  Figure  18.  From  the  diagram. 


V 


observed 


where 

RISW  =  RIS  (wall) 

RISP  =  RIS  <Profile) 

RISp  =  fj  (emf) 

RISW  =  f2  <emf) 

From  the  observations  reported  in  the  previous  sections,  it  can  be  concluded  that  f2  is  a 
stronger  function  than  fj.  since  an  increase  in  wall  leakage  resistance  seems  to  accompany 
the  switch  from  the  nonsaturated  to  the  saturated  mode. 

Before  a  final  analysis  can  be  made,  a  more  extensive  study  must  be  undertaken  and  additional 
data  collected  regarding  leakage  behavior  so  that  the  functions  fj  and  f2  can  be  firmly  established. 
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VI.  NON  EQUILIBRIUM  PROCESSES  IN  MAGNE  I  OACTIVE  PLASMA  E  LOWS 

hi  a  plasma  to  which  an  electric  field  is  applied  approximate  local  statistical  equilibrium  can 
be  assumed  for  each  component  species  when  the  energy  acquired  by  an  elei  tron  along  a  mean 
free  path  is  small  — by  a  factor  of  order  (n,eieetroii/  rnatom^ 1 '  compared  to  the  thermal 
energy.®  Howevir,  the  thermal  equilibrium  between  the  electron  component  and  the  atomic 
components  is  frequently  disturbed,  as  the  exchange  of  energy  in  elastic  collisions  between 
the  light  electrons  and  heavy  atomic  particles  is  extremely  small.  The  nonisothermal  plasma 
state  is  generally  associated  with  an  unusually  high  electron  density  indicating  a  quasi-ioni¬ 
zation  equilibrium  with  the  electron  temperature.®  This  is  observai  ie  in  low  pressure  dis¬ 
charges  (m.  f.  p.  »  plasma  dimension,  recombination  at  wall)9  as  well  as  in  hig..  pressure 
discharges  (m.  f.  p,  «  plasma  dimension,  recombination  in  volume)10 

Instead  of  applying  an  electric  field  to  a  plasma  at  rest,  an  elect,  ic  field  action  can  be  induced 
by  moving  the  plasma  across  a  region  occupied  bv  a  magnetic  field,  where  by  means  uf  an 
appropriate  electrode-circuit  system  the  closing  of  the  induced  currents  is  provided.  This 
leads  to  the  concept  cf  magnetoactive  plasma  flow.  In  contrast  to  the  resting  systc-m  in  which 
unlimited  time  is  available  for  the  development  of  the  nonequilibriuni  state,  the  time  available 
in  the  steady  flowing  system  for  the  development  of  any  nonequilibrium  processes  is  limited  to 
the  interval  required  for  a  plasma  volume  element  to  cross  the  magnetic  field  region.  The 
plasma  in  the  volumg^element  changes  its  properties  with  position  aLo  due  to  the  effect  of  the 
magnetic  body  force  and  the  gasdynamir  forces,  Accordingly,  an  observer  moving  with  an 
arbi  rary  plasma  volume  element  would  observe  that  the  fields  change  with  time.  Considered 
from  the  moving  system  of  reference,  a  steady  state  is  not  reached  before  the  plasma  volume 
element  has  left  the  magnetoactive  region  and  the  fields  have  relaxed  to  their  thermodynamic 
equilibrium  values. 

In  the  following  an  attempt  is  made  to  investigate  theoretically  the  complex  and  interacting 
processes  in  subsonic  plasma  flows  transverse  to  homogeneous  exterior  magnetic  fields.  The 
main  attention  is  given  to  the  processes  which  lead  to  the  buildup  of  nonequilibrium  electron 
temperature  and  nonequilibrium  electron  density.  The  analvsis  is  based  on  multifluid  field 
equations  derived  from  the  Boltzman  equation  in  the  13-moment  approximation. 

FIELD  EQUATIONS 

The  following  analyses  consider  steady  flows  (d/dt  =10)  cf  low  temperature  plasmas  consisting 
of  a  driving  gas  (o)  with  high  ioniz,  ion  potential  (in  which,  consequently,  reactions  can  be 
neglected)  and  the  ionization  products  of  an  easily  ionizable  seeding  gas  — seed  atoms  (a),  seed 
ions  (1),  and  electrons  (e).  In  particular,  the  numerical  analysis  refers  to  a  1%  cesium  -  seeded 
helium  plasma  initially  at  a  pressure  of  PQ  =  1.5  X  10  new  ton /m*5  and  a  temperature  of  T0  = 
200()°K  (o  =Hr,  a  SCs,  i~Cs+).  For  electron  temperatures  of  the  same  order  of  magnitude 
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as  the  gas  temperature,  ionizations  and  recombinations  in  stops  can  be  disregarded  in  the  first 


approximation.11  For  the  plasma,  the  essential  reaction  mechanisms  are  two -body  electron- 
atom  collision  ionization1* and  three-body  cleetron-ion-electron  collision  recombination,13 
from  and  to  the  atomic  ground  level,  respectively.  The  corresponding  reaction  equation  \s: 


e  +  at^i  +  e  +  e 


(1) 


Introducing  the  stoichiometric  coefficient,  vs,  of  the  components  s  =  e,  i,  a,  and  o  in  accordance 
with  Equation  (l) 


ve  =  +  2  -  1  =  +  1, 


l  T 


•'o'  0 


(2) 


the  reactive  mass  density  production  in  the  components  (s)  is  described  by 


u 


*  Tns  ms  +  ns  mfiv  *  v*c  =  vs  ins  (  -  I’  ) 


(3) 


The  reaction  rates  for  ionization,  1’  1 ,  and  rc  combination,  r~,  are  related  to  the  particle  den¬ 
sities  and  electron  temperature  by 


_  n  Q 

1  ne  na  ea* 


f  np  rij  np  Reje 


(4) 


The  probability  coefficients  for  two-body  (e-a)  ionization,  Sea,  and  three-body  (e-i-e)  recombi¬ 
nation,  R  •  ,  for  Maxwellian  plasmas  are:15* 
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(6) 


where  £  =  number  of  electrons  in  the  highest  shell  n  of  the  atom  in  ground  state,  =  ioni- 

zation  energv  of  the  atom  from  n-th  shell,  U.  ;  =  2  J  °'.  +  1  =  excitational  partition  functions. 


and  G  = 


*For  small  seeding  ratios,  (nj  +  na)  <nQ,  the  perturbation  of  the  electron  velocity  distribution 
by  the  reactions  can  be  disregar  led.1* 
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As  a  conclusion  from  Equation  (3),  the  total  mass  density  is  conserved  in  the  reactions, 

(I'+  -I’  )  0.  The  total  particle  density  is  not  conserve  d  in  the  reactions, 

s 

y>s  ( T+  -  1"')  4  0,  the  trivial  case  of  reactive  equilibrium,  T  '  T  ”,  being  excepted. 


When  the  plasma  enters  the  region  occupied  by  the  transverse  magnetic  field,  the  charged 
particles  are  diverted  into  a  gyrational  motion  around  ttie  magnetic  field  lines,  while  the  neutral 
particles  are  influenced  indirectly  by  the  magnetic  field  via  collisions  with  the  charged  particles. 
As  a  consequence  of  the  collisions  of  the  charged  particles  with  the  neutral  particles  and  with 
themselves,  tneir  gyration  centers  are  displaced  in  a  well  defined  way  parallel  and  transverse 
to  the  flow.  In  the  macroscopic  picture,  this  means  a  drift  of  the  plasma  components  relative 
to  each  other  builds  up  spatially  in  these  directions.  For  the  component  (s),  whose  particle 
interactions  with  the  remaining  components  (r  4  s)  result  macroscopic-ally  in  intercomponent 
friction  forces  [~(v^  -  vp],  this  process  is  described  by17 


nsms  vs 


V  vB 


VPs 


(E  + 


vs  X 


B) 


n 


s  ms 


E 

r  i  s 


-1 

sr 


(Vn  -  Vr) 


(7) 


The  various  combinations  of  interactions  between  the  charged  and  neutral  particles  in  plasma 
are  treated  as  Coulomb  interactions  I  Debye-radius  D  =  (4w^ns  eg  /4w  <  0  k  Ts)”  as  maxi- 

L  n  s  2 

mum  impact  parameter;  transport  cross  section  QS1.  =  —  (es  er/4nt0kT  )4  •  lnAsr; 

7|os  er,  /  \-»1"  "  ‘  f 

Asr  -  Dl -  /  .I  kT srj  and  elastic  sphere  interactions  [interaction  radii  rs  and  rr; 

nt°  '  '  o  ,»• 

transport  cross  section  Qsr  =  tr  frs  +  rr)  J  respectively.  The  times,  t  sr,  determine  the 
relaxation  of  the  velocity  fields,  and  are  given  by 


-l  =  1J1 
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mgr 
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nr  Q 


r  ^sr 


(8) 


where 

msr  *  msmr/(ms  +  mr)  =  reduced  mass 

Tsr  =  msr  [(Ts/ms)  +  (Tr/mr)]  =  reduced  temperature 

The  electromagnetic  field  B,  E  acting  on  the  plasma  is  composed  of  the  exterior  fields  (sub¬ 
script  "o")  and  the  fields  produced  by  the  motions  of  the  charges  (subscript  "s"): 

B  =  B(J  *■  B,|t,  etc.  For  a  quasi-neutral  plasma,  E  =  E  ,  and 

vxe  |*0J]nsPsvs,  V-  1  =  0  (P) 

s  ° 
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V  X  E  =  fi,  T  •  P!  =  0  (10) 

The  omission  of  space  charge  effects  is  based  on  the  presumption  that  the  characteristic  spatial 
extension  of  the  plasma  is  large  compared  to  the  Debye  shielding  distance,  D. 

For  the  subsonic  low  temperature  flows  under  consideration,  energy  transport  by  radiation  and 
thermal  heat  fluxes  [q^  -  q*st  a'*ia^ve  t  ^thermal  ]  is  not  of  primary  importance.  *  The 
significant  irreversible  processes  are  intercomponent  frictional  heating  [<~(v^  -  v*r)^]  ,  inter¬ 
component  scale.  heat  flux  [~(TS  -  T  J],  and  reactive  energy  production  [~(T+  -  F-)].  Under 
formal  inclusion  of  the  sources  of  qs,  the  change  of  energy  density  in  the  component  (s)  is 
described  by17 


[%•  VPS  +7S  Ps  v-  vs]  (Ts  -  1) 


-1  V-'  Ts  msr  -l  -*  >2 

a  ~ - —  T  sr  *  ns  ms  <vs  -  V  " 


r  /  .  T« 


3  Z 


m 


sr  - 1 


4  s 


m. 


r 

sr 


ns  k  (Ts  -  T ,,)  -  *s  *  %  <r*  -  r‘)  -  V  •  qg 


(11) 


In  an  ionization  act,  the  ionization  energv  is  taken  out  of  the  electron  component.  In  a  recom¬ 
bination  act,  the  binding  energy  of  the  same  amount  is  fed  into  the  electron  component.  Fur¬ 
thermore,  with  every  transformed  s-particle,  the  energy  3/2  k  Ts  is  liberated  in  the  mean. 
Therefore,  in  Equation  (11),  the  reactive  variation  of  energy  per  particle  of  the  component  (s) 
is  given  by  (8rs  =  Kronecker  8  -  svmbol) 

^  *  s  ~  *  n  ® es  ”  ^ 


By  summation  of  Equation  (It)  over  s  it  can  fie  shown  —  with  help  of  Equations  (3)  and  (12)  —  that 
the  plasma  as  a  whole  exhibits  a  reactive  energy  source,  S  =  «  n  ( r4  -  F  ). 

The  equations  of  state  for  the  individual  components  (s)  are: 

pg  =  ns  k  Tg  (13) 

Equatic  s  (3),  (7),  (9),  (11),  and  (13)  form  a  complete  set  of  field  equations  for  the  considered 
steady  multitemperature  multicomponent  flow.  The  expressions  for  the  elastic  transport  pro¬ 
cesses  are  valid  for  infrasonic  relative  drift  velocities.17 


1  _  2 

msr  (vs  -  vr)  «  k  Tsr  (14) 


-  A  discussion  cf  the  assumptions  made  is  presented  later,  together  with  furthei  problems. 
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PLASMA  FLOW  MODEL 


Consideration  is  given  to  a  parallel  plasma  flow  of  initial  velocity,  VQ,  (I  (-direction)  entering 
at  x  =  0,  a  region  occupied  hv  a  transverse  homogeneous  induction.  B0,  extended  up  to  x  L. 
In  the  directions  perpend:  uiar  to  V0  and  B0  (l-direction),  the  flow  is  limited  by  electrodes 
[parallel  to  the  plane  (VQ,  Bn)  and  a  distance,  d,  apart]  which  collect  the  induced  electrical 
currents  into  an  exterior  circuit  system.  In  the  directions  parallel  to  BQ,  the  flow  is  assumed 
to  be  unlimited. 

For  mathematical  convenience,  the  following  conditions  are  introduced,  allowing  the  flow  to  be 
treated  as  quasi-one-dimensional  [in  which  any  flow  field  can  be  only  a  function  of  the  axial 
coordinate,  x:  fs  =  fg  (x)].  The  magnetic  stresses  are  assumed  to  be  large  compared  to  the 
viscous  stresses.  This  requires,  or  V  B^»ftV  d*^  —  i.e.,  Hartmann  numbers20 


H  =  (M/tr)1  /2  d  B  »  1 


(15) 


Under  this  condition,  the  transverse  inhomogeneitios  due  to  viscous  stresses  can  be  disregarded 
and  the  nonuniformity  in  the  axial  direction  is  influenced  solely  by  the  force  density,  k(|  =  j  ^ 

X  B,  resulting  from  the  interaction  of  the  transverse  current  density,  ^  ns  es  VsJL>  with 

s 

B.  At  inductions  B0,  for  which  w  e  r  p  >1,  o>c  -  |ee|  B/me,  and  re“*  =  ^  res»  the  resul~ 

s  i  e 

tant  of  the  Hall-field  and  axial  electron  pressure  gradient  diffusion  field  would  produce  a  down¬ 


stream  current  dens 


i.y,T„  =£> 


_  e  v  .  The  interaction  of  the  latter  with  B  would  result  in 
S  s  s  U 


a  transverse  force  density,  k  =  j(|  X  B,  perturbing  the  one-dimensionality.  To  ascertain  the 
one-dimensionality,  the  resultant  axial  plasma  field  is  assumed  to  be  compensated  by  r  field 


{j„ 5  «■ 


[E-s 


e  ~v  .  X  B  -  v  p  I  (n„  |e  I ) 
s  si  v ir  e  I  e |  el 


i-l 


(16) 


This  requirf  an  ideally  segmented  electrode  system,  each  opposite  electrode  pair  being  con- 

21 

nected  bv  an  individual  exterior  circuit.  Because  of  Equation  (16),  the  electrical  power  liber- 

ated  in  the  plasma  is  that  of  an  ideal  Faraday  flow,  W  =  j  •  (E  +  V  X  B)  =  <r(E  +  V  X  B) .  , 

2 

where  <r  =  ne  ee  r  e/me  is  the  scalar  conductivity.  So-called  end  losses  and  their  effects  are 
ignored.  These  can  be  estimated  separately .** 


For  an  analysis  of  nonequilibrium  effects,  it  is  adequate  to  allow  for  a  maximum  induced  e!e<  - 
trie  field  action.  For  this  reason,  opposite  electrodes  are  taken  as  short  circuited  by  their 
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exterior  circuit.*  Th«  geometry  of  the  flow  fields  can  now  be  stated  at-  foil  ws  for  the  eon 
side  red  one-dimensional  flow  model  (Figure  10): 


h  ° 


R  =  B  +  R  ; 
o  * 


B*  I!  B0 


vs  vS)|  +  vSJ_, 


ve  1 1  vi  i 


-x  =  0  x  =  x 


x  >0 


x  =  L 


4400-16 


Figure  19.  Geometry  of  flow  fields. 


For  the  numerical  analysis,  a  plasma  is  specified  which  is  (  imposed  of  helium,  cesium,  and 

-2 

the  ionization  products  of  the  latter  at  a  se^d  ratio  a  =  (na  n^)nQ  =  10  .  The  initial  ioni¬ 

zation  is  taken  in  accordance  with  the  Saha  equation  (thermodynamic  equilibrium).  The  initial 
values  defining  the  flow’  before  entering  the  magnetic  field  region  are  assumed  as  shov/n  in 
Table  II.  Different  exterior  inductions  are  considered  — 130  =  2,  3,  4,  and  5  Tesla. 

Table  II. 

Initial  flow  values. 


Pressure  —  PG  =  1.5  X  103  newton/m^ 
Mach  number — M0  =  0.5 


Temperature  —  Tt)  =  2000°K 
Velocity— VQ  -  1.155  X  103  m/sec 


Helium  density — nQ  0  =  5.38  X  103^  m'3  Cesium  density — na^  {)  =  5,38  X  1033  m  3 

Electron  density — n  =  4.24  X  10^  m-3  Ion  density — n-  =  4.24  X  10*°  m-3 
J  e,  o  l,  o 


24  1 

The  elementary  constants  of  the  plasma  particles  are  listed  in  Table  III  (Qsr  =  Qrs).  ' 


*The  case  of  ohmic  exterior  circuits,  w'hich  is  primarily  of  technical  interest,  has  been  studied 
in  an  earlier  investigation. ,s 
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Table  III. 

Elementary  particle  constants. 

Qpo  =  5.5  X  10"20  m2  mQ  -  6.540  X  10'27  kg 

Qoa  4.0  >  10"18ni2  nia  2.172  X  I0‘25  ki> 

Qio  9,0  X  10'20m2  uij  2.172  X  tO-25  kg 

Qia  2.0  X  10" 18  m2  ni0  =  9.108  X  10"31  kg 

«n  3.893  eV  e  =  1.602  X  10"19ampsec 

No  evaluation  is  made  of  the  total  generated  current  (current  density  integral  over  the  electrode 
plane).  In  this  case,  the  n  meric al  analysis  does  not  require  specification  of  the  electrode 
dimensions,  nor  of  the  electrode  distance  as  the  exterior  electric  field  is  given. 

NUMERICAL  RESULTS  AND  ANALYTICAL  DISCUSSIONS 


The  field  equations  of  the  one-dimensional  flow  model  were  integrated  numerically  for  the  given 
initial  conditions.  The  computer  program  employed  a  second  order  Runge-Kutta  numerical 
integration  scheme  to  carry  the  investigation  over  the  entire  subsonic  flow  length  —  i.e.,  up  to 
that  point  where  the  Mach  number  becomes  M  =  1.  In  the  regions  of  strong  field  gradients  an 

_  c 

integration  step  size  of  Ax  -  10  cm  was  used  to  ensure  uniformity  and  accuracy. 

Change  of  Mean  Flow  Fields 


The  interaction  of  the  induced  jurrent  density  with  the  magnetic  induction  results  in  a  force 
densitv  antiparallel  to  the  flow  velocity,  k^  =  -  B  e)(,  Therefore,  the  magnetic  interaction 
has  an  effect  on  the  flow,  which  is  quantitatively  simila  to  a  viscous  drag  force.  This  effect 
leads  to  the  terms  "magnetic  drag  force”  and  "magnetic  viscosity.  "  Using  this  analogy,  it  is 
concluded  that  the  mean  plasma  fields  change  with  increasing  length  coordinate,  x  >  0,  in  a 
subsonic  flow  as  follows: 


«)  Mean  static  pressure,  P  =  ps,  mean  static  temperature  "  =  ns  7s^Z»ns'  and 

s  s  s 

mean  mass  density,  p  -  nsms,  decrease. 


/S)  Mean  mass  velocity,  V  -  ns  ms  v^/^  ns  ms,  mean  dynamic  pressure,  Pd  =  —  n 

s  s  2  s 

and  Mach  number,  M  =  |  ns  msvs  j/^^^s  ng  k  Ts,  increase. 


"  msvs  , 

&  s 


Accordingly,  under  subsonic  conditions,  ultimately  a  transformation  of  potential  energy  or 
static  pressure  into  kinetic  energy  or  dynamic  pressure  takes  place — in  the  case  of  an  ohmic 


37 


.AlliMon 


exterior  circuit,  electrical  energy  would  also  he  gained.  In  Figure  20,  the  quantitative  develop 
ment  of  the  dynamic  processes  a  -yS  is  shown  in  the  considered  flow  model,  where  B()  is  a 
parameter.  Note  that  the  length  at  which  the  flow  becomes  soni  M  1,  decreases  strongly 
with  increasing  induction,  BQ. 


The  flow  behavior  can  be  understood  physically  from  the  second  law  of  thermodynamics.  For 
a  simple  explanation,  eliminate  from  the  field  equations  in  the  one-fluid-approximation  the 
gradient  of  the  flow  velocity  (*  =  E  /VB  =  load  factor):* 


(20) 


where  for  an  adiabatic  flow,  1  =  0,  and  for  a  nonadiabatic  flow 


1  = 


<y  -  i) 


■> 

(i  -*)“ 


*  n<r  -  r~)  +  v 

9  9 

if  V"  B1" 


1  0 


(21) 


For  the  flow  under  consideration,  the  factor  [(1  -  k  )  +  is  necessarily  positive.  It  follows 

M  <  1:  cIV/dx  >0;  M  >  1:  d  V./dx  <  0 

Whether  a  ansformation  of  potential  energy  or  kinetic  energy  takes  place  in  the  magnetoactive 
region  depends  consequently  on  whether  the  flow  is  initially  subsonic  or  supersonic. 


^Summation  of  Equations  (7)  and  (11)  over  s  =  e,  i,  a,  and  o  yields: 

PV  •  VV  =  -  VP  +  7*  B  +  A,  . 

v  ■  v  p  +  yp  v  •  v  -  ( r  -  1 ) [V *  (e  +  v  x  b)  -  <  n  (r‘  -  r')  -  v  •  "qjt  a2 

It  is  usual  to  neglect  the  small  quantities  A^,  and  A^  (wg  =  vg  -  V): 

41  ’E  "sms''s  • 


-E 


S  c 


w 


TPs  +rKPs  v  ws  -—  (ys  -  1)  rs  kTs  (r+  -  r') 


Equation  (20)  follows  then  from  (A)  and  (B)  by  elimination,  using  Ohm’s  law  in  addition. 


(A) 

(B) 
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Figure  20.  Mean  flow  fields  (upper  abscissa  for  BQ  =  2  Tesla) 


.  Allison 


Integration  of  Equation  (20)  from  x  =  0  (V  -  VQ,  Tvl  M0)  to  x  =  L  (V  =  Vj,  M  -  Mj)  loads  to 
the  following  integral  equation  for  the  length,  L,  at  which  Hie  Mach  number  becomes  Mj,  where 
consideration  is  made  of  the  invariance  of  p\ 
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o  v  o 
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f  <t  [(1  -*)  +  S]  B2  dx 
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(22) 


Adapting  to  the  conditions  of  the  flow  model  by  putting  *  =  0, 
the  left  side  of  Equation  (22)  in  quasi-sdiabatic  approximation 
where  y  *  is  a  quasi-adiabatic  coefficient]  ,  yields: 
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B  -  -  1),  and  integrating 

[d  VIV  =  2[y  */{y*  +  1)]  dM/M, 


(23) 


where  a  is  a  mean  conductivity  defined  bv: 


I 


0 


dx 


(24) 


For  a  given  set  M  .,  L  is  inversely  proportional  to  the  mean  conductivity,  ,  and  inversely 
proportional  to  the  square  of  the  induction,  BCj .  It  must  be  noted,  however,  that  cr  depends 
indirectly  on  BQ,  These  elementary  considerations  are  reflected  in  the  numerical  results. 
Figure  20, 


Buildup  of  Nonequilibrium  Electron  Temperature 


The  directed  kinetic  energy  acquired  by  the  charged  particles  in  the  induced  electric  field  is 
transformed  in  part  into  heat  by  intercomponent  friction.  As  the  fraction  of  relative  energy 
lost  by  the  s-particles  in  elastic  collision  with  r-particles  is  of  the  order  ms/(ms  +  mt>),  the 
thermal  equilibrium  between  the  various  heavy  particles  is  only  weakly  disturbed.  The  elec¬ 
trons,  however,  can  dispose  of  the  steady  heat  input  only  if  their  energy  is  relatively  great 
compared  to  that  of  the  heavy  particles  as  the  fraction  of  relative  energy  lost  in  elastic  collisions 
with  heavy  r-particles  is  of  the  order  me/(me  t  mr).  Therefore,  the  electron  gas  heats  up 
relative  to  the  remaining  plasma  until  the  losses  in  elastic  collision  with  particles  of  the  other 
components  r — which  result  macroscopically  in  scalar  heat  flows  ~(Te  -  Tr)  —  and  other  loss 
mechanism  ,  balance  the  electron  heating.  This  balance  is  established  at  a  length,  !  j.,  which 
does  depend  explicitly  on  the  induction,  Bn.  Above  this  length,  x  >t  -p,  the  flow  exhibits 
a  practically  constant  electron  temperature  plateau  up  to  that  point  where  the  Mach  number 
begins  to  rise  steeply,  M  -*l.  In  this  region,  due  to  the  strong  increase  of  induced  electric 
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field  arid  electron  heating,  the  electron  tcmnerature  lumps  abruptly.  Figure  21  shows  quan¬ 
titatively  how  the  electron  temperature  builds  up  in  the  considered  flow  model  at  the  various 
inductions,  B0.  Note  that  (T{5  -  Tn)  increases  at  a  rate  less  than  proportional  to  B^  lor  increas- 
mg  B0. 


The  buildup  of  the  intercomponent  thermal  nonequilibrium  can  be  explained  on  an  elementary 
analytical  basis  by  considering  that  within  0  <  x  £  t.p  tin  fields  —  with  exceptiorr  of  Tp —  remain 
practically  constant,  and  that  the  electron -neutral  interactions  predominate  under  the  c.ondi- 
tiona  given.  Disregarding  the  strong  T0 -dependence  of  the  Coulomb  interactions,  (Qpl  ~  Tp  ) 
Equation  (11),  taken  for  s  -  e,  indicates  a  convective  buildup  length  of  the  order 
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(25) 


Neglecting  the  Coulomb  interaction  term,  one  finds  with  the  initial  values  and  constants  of 
Tables  II  and  III 


Above  x  =  j6^,,  the  plasma  is  quasi-homogeneous.  Consequently,  in  establishing  the  energy 
balance  for  the  electron  component  in  the  plateau  region,  the  convective  terms  in  Equation  (11), 
taken  for  s  =  e,  may  be  neglected.  The  balance  gives  —  after  eliminating  in  the  same  approxi¬ 
mation  the  drift  velocities  by  means  of  Equation  (7)  —  for  the  plate  au  electron  temperature 

<Ti,a  =  To> 


where  +  V  X  B  is  the  field  producing  the  electron  heating  and  E^  is  the  critical  plasma 

field  for  electron  heating,7 
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The  last  term,  8  ,  is  due  to  the  inelastic  energy  source,  and  of  importance  in  strong  reactive 
nonequilibrium,  P  +  i  F”, 
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(29) 


This  equation  is  transcendent  in  T„,  and  admi:s  of  no  exp’icit  analvtical  solution.  However, 
neglecting  the  Coulomb  interaction  terms  and  the  inelastic  source  terms.  Equation  (20)  becomes 
a  quadratic  equation  in  Te/T  .  The  latter  has  the  solution: 


1  +\/l  +  4  (E'jEfy 


(301 


where 


E.  --  — , 


1 - 

r  "I 

-  ** 

1  2  me 

j  i  mo  na  ^ca 

.  n  a  ^ea 

1 3«r  m0 

nia  no  Qeo 

no  ^eo 

' 

L  J 

L-  -J 

no  ^  T0  *  Qoo 


The  initial  values  and  constants  of  Tables  II  and  111  give 


Ej  =  29.31  v/cm. 


(31) 


As  an  application,  the  plateau  electron  temperatures  are  given  in  Table  IV  according  to  Equation 
(30)  for  the  different  inductions  BQ. 


Table  IV. 

Plateau  electron  temperatures. 


Bo- 

Tesla 

2.00 

3.00 

4.00 

5.00 

A- 

V  Br. — v/cm 
o  ° 

23.10 

34.65 

46.20 

57.75 

Te 

°K 

2866 

3565 

4309 

5070 

The  deviations  from  the  numerical  Te-values  (1  to  10%)  indicate  tha.  at  extremely  high  induc¬ 
tions  thf»  effects  of  flow  inhomogeneity  and  reactive  nonequilibrium  acquire  increasing  influence. 

Buildup  of  Nonequilibrium  Eleetron  Density 


To  understand  the  ionization  in  a  steady  plasma  flow,  consider  first  a  highly  idealized  map- 
netoact  ve  plasma  flo'A.  The  latter  is  quasi-homogeneous  and  may  exhibit  a  relatively  infinitely 


43 


\ili«*on 


extended  electron  temperature  plateau,  T  =  constant  >  T0  above  a  negligible  small  length, 

The  electron  density  at  x  *  0  is  in  equilibrium  with  the  gas  temperature  np  0  n(1  ()  (TQ). 

Further  aSaunn  that  the  develop:-  ionization  reactions  do  not  disturb  the  uniformity  >f  the  flow 
and  that  the  electron  density  remains  small  com  oared  to  the  atom  density  n<.«  na»  i.e* 
nQ  Sna  0  constant.  From  Equation  (3),  taken  for  s  =  e,  there  follows  under  these  conditions 
the  "relaxation"  length  within  which  the  electron  density  builds  up  due  to  ionization: 


*8 
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na,  o  ,\‘a  ^  *  e^ 


132) 


In  Table  V,  (  s  is  given  based  on  the  initial  values  and  constants  of  Tables  II  and  Ill,  while  Tp 
is  taken  in  accordance  with  the  plateau  electron  temperature  of  Figure  21. 


Table  V. 

Relaxation  lengths. 


B0  —  Tesla 

2. 00  3.0.:  4.00  5.00 

TC~'K 

2820  3455  4070  4615 

( s—  cm 

2. 73  X  102  i . 36  X  1 0 1  1. 76  X  10°4. 45 X  10' 

In  Figure  22,  it  is  shown  quantitatively  how  the  nonequilibt  ium  electron  density,  ne,  develops 
in  the  considered  flow  model  at  different  inductions  B0.  As  an  essential  result,  we  recognize 
that  ne  does  not  increase  notic  eably  in  the  Te-plateuu  region.  For  explanation,  compare  the 
hypothetical  ionization  length  £  ,  given  previously,  with  the  extension  of  the-  Tc-plateau,  to  be 
estimated  from  Figure  21.  Throughout  the  Tp-plateau  region,  nr  does  not  correspond  to  an 
ionization  equilibrium  determined  by  Te,  rather  in  orders  of  magnitude  0[np]?0  [nn  G] . 

Under  the  strong  nonequilibrium  conditions  developing  at  the  pnd  of  the  subsonic  flow  zone, 
r+  and  n(,  increase  a  ten  ply,  however,  even  there  no  ionization  equilibrium  determined  by  Tp 
is  reached.  For  explanation,  it  is  remembered  that  the  change  of  Te  is  produced  bv  elastic 
collisions  and  accordingly  represents  a  relatively  quick  process,  while  the  change  of  ne  is  pro- 
d  iced  in  the  main  by  inelastic  collisions  and  accordingly  represents  a  relatively  slow  process. 

An  increase  of  the  electron  dens.ty  would  increase  the  magnetic  force  density  and  make  the  flow 
jump  to  M  — »  1.  It  appears  that  these  effects  —  increasing  nf,  and  M  —  condition  each  other. 

This  would  mean  that  the  result  obtained  for  our  special  flow  model  —  namely,  that  electron 
density  increases  strongly  in  the  final  region  where  the  flow  becomes  sonic,  has  a  much  more 
general  validity. 
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Finally,  in  Figure  23,  it  is  shown  quantitatively  how  the  current  density,  j,  changes  in  the  con¬ 
sidered  flow  model  for  different  inductions,  RQ.  It  is  noteworthy  that  in  the  region  of  the  elec¬ 
tron  temperature  plateau,  j  is  depressed.  This  is  due  to  the  reduction  of  the  electrical  con- 
ductivity  by  the  elevation  of  Tc  at  constant  ne  (c  ~  ne  Te  Beyond  the  end  of  the  Tp-plateau 

region,  the  electrical  conductivity  increases  due  to  the  stronger  dependence  on  ne.  As  the 
induced  electric  field  increases  in  this  region,  j  increases  there  also.  The  magnetic  force 
density  changes  proportional  to  j,  as  the  induced  magnetic  field,  B^,  is  small  compared  to  the 
applied  field,  BQ. 

CONCLUSIONS 

A  model  for  a  one-dimensional  magnetoactive  flow  of  a  nonisothermal  reacting  four-component- 
plasma  has  been  conceived.  Based  on  multifluid  field  equations,  the  nonequilibrium  processes 
developing  in  the  subsonic  region  of  the  flow  have  been  computed  via  numerical  integration. 
Parallel  to  that,  an  elementary  analytical  theory  has  been  established  explaining  the  numerical 
results  in  terms  of  simple  functional  relationships.** 

Understanding  that  quantitative  statements  apply  to  the  considered  plasma  model  and  initial 
conditions,  the  following  conclusions  can  be  drawn. 

•  The  force  density  r  esulting  from  the  interaction  of  the  induced  current  density  and  the 
magnetic  field,  and  the  nonadiabatic  effects  cause  a  decrease  of  the  static  pressure  and  an 
increase  of  the  dynamic  pressure  in  the  subsonic  flow.  These  field  changes  are  most  pro¬ 
nounced  close  to  the  end  of  the  subsonic  flow  region  where  the  singularity  M  =  1  is  approached, 

•  The  power  liberated  by  the  friction  of  the  interpenetrating  components  causes  the  flow  to 
becor.  s  nonisothermal,  the  nonisothermv  being  quantitatively  important  between  the  elec- 
tron  gas  and  the  remaining  plasma.  After  a  buildup  length  of  f -p  =  10  cm,  the  electron 
component  exhibits  an  elevated  nearly  constant  temperature  (Te-plateau)  extending  up  to 
that  point  beyond  which  the  Mach  number  increases  strongly  towards  M=l. 

•  The  electron  density  does  not  change  its  order  of  magnitude  (compared  to  its  initial  value) 
throughout  the  Tp-plateau  region,  even  at  extreme  high  BQ-fields,  The  flow  exhibits  a 
strong  reactive  nonequilibrium  in  which  the  electron  density  does  not  correspond  to  its 
equilibrium  value  with  the  electron  temperature.  The  assumption  of  an  ionization  equili¬ 
brium  in  that  region  made  frequently  in  the  literature  must  be  rejected. 

•  Close  to  the  end  of  the  subsonic  flow  region  the  Mach  number  increases  steeply  towards 
the  singularity  M=l,  where  the  nonequilibrium  conditions  become  extremely  strong.  Elec¬ 
tron  temperature,  ionization  reaction  velocity,  and  electron  density  practically  jump  in 
this  region.  There  the  electron  density  changes  its  order  of  magnitude.  Recombination 
processes  are  unessential  in  the  whole  subsonic  flow  region  (ne  is  always  below  its  equili 
brium  value  with  Te,. 
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The  calculations  have  been  carried  through  in  a  dissipative  approximation  considering  inter¬ 
component  fric  lon.  Consequently,  the  extremely  strong  nonequilibrium  effects  obtained  close 
to  the  end  of  the  subsonic  flow  < M  — *  1 )  must  be  interpreted  more  qualitatively  than  quantitatively. 
Due  to  the  strong  inhomogeneity  of  the  flow  in  this  final  region,  viscous  dissipation  (~V“  vs)  and 
thermal  conductive  dissipation  Ts)  acquire  increasing  influence,  and  have  to  be  taken  into 
account  in  a  more  refined  theory. 
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VII.  CONCLUSIONS 

During  the  report  period,  a  closed  loop  MPD  power  generation  device  was  designed,  constructed, 
and  operated.  Two  modes  of  operation  were  observed  —  saturating  and  nonsaturating.  The 
saturating  mode  is  associated  with  low  seeding  ratios  while  the  nonsaturating  mode  is  associated 
with  high  seeding  ratios  and  high  leakage  resistances.  During  operation,  leakage  resistances 
varied  between  0.1  and  1.0  megohms. 

Of  the  two  modes  observed,  only  the  saturating  mode  was  investigated  thoroughly.  The  obser¬ 
vations  were  made  under  the  following  conditions: 

•  Pressure — slightly  above  atmospheric 

•  Magnetic  field  —  up  to  21,  500  gauss 

•  Velocities — up  to  250  m/sec  at  the  channel  entrance 

•  Temperatures  (static)— up  to  1500TC 

For  the  saturating  mode  at  these  conditions,  no  nonequilibrium  ionization  was  observed. 

Theoretical  considerations,  as  discussed  in  Section  VI,  indicate  that  to  observe  nonequilibrium 
ionization,  runs  should  be  made  at  lower  pressures  and  higher  velocities  and  temperatures 
than  those  stated  previously.  Such  runs  are  planned  during  the  remainder  of  the  program. 
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